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Abstract : This paper deals with homogenization of second order divergence form 
parabolic operators with locally stationary coefficients. Roughly speaking, locally 
stationary coefficients have two evolution scales: both an almost constant micro- 
scopic one and a smoothly varying macroscopic one. The homogenization proce- 
dure aims to give a macroscopic approximation that takes into account the micro- 
scopic heterogeneities. This paper follows |[T3ll and improves this latter work by 
considering possibly degenerate diffusion matrices. 

Resume : Nous etudions l'homogeneisation d'operateurs paraboliques du second 
ordre sous forme divergence a coefficients localement stationnaires. Ces coeffi- 
cients presentent deux echelles d' evolution: une evolution microscopique presque 
constante et une evolution macroscopique reguliere. Latheorie de l'homogeneisation 
consiste a donner une approximation macroscopique de l'operateur initial qui tient 
compte des heterogeneites microscopiques. Cet article fait suite a |[T3l et generalise 
ce dernier en considerant des matrices de diffusion pouvant degenerer. 

AMS classification: 60F17; (35B27; 35K65; 28D05). 



1 Introduction 

This paper follows fT3l and deals with homogenization of second order PDEs with locally sta- 
tionary coefficients by means of probabilistic tools. More precisely, we aim at describing the 
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asymptotic behavior, as e goes to 0, of the following Stochastic Differential Equation (SDE) 

rt ve i-t ve i-t 



(i) x, 



- f b(u^,X £ r )dr+ [ c (uj,^,X £ r )dr+ [ a(u, ^, X £ r ) dB r , 
e Jo £ Jo £ Jo £ 



where B is a standard d-dimensional Brownian motion and the parameter uj evolves in a random 
medium 0, that is a probability space with suitable stationarity and ergodicity properties. For 
each fixed value of the parameter y G R. d , the coefficients b(u>, -,y), c(u>, -,y) and a(uj, -,y) are 
stationary random fields (the parameter to stands for this randomness). That is why they are said 
to be locally stationary. The generator C e of the process X s can be written in divergence form 
as 

i,j=i j 

for an antisymmetric matrix H, a real- valued function V and a = aa*. 

Let us first briefly outline the chronological approach of this issue. The convergence of the 
previous SDE (or the connected PDE) has been first established in the locally periodic case, that 
is when the coefficients are deterministic and periodic with respect to the variable x/e (l]l2l- 
Due to the lack of compactness of a random medium, the random case raises more difficulties. 
As far as we know, the first work in this context is due to Olla and Siri in ifTTTl . The authors 
considered a nearest neighbors random walk on 7L evolving in a locally stationary environment. 
They established an invariance principle for this process under diffusive scaling of space and 
time. The main tool of the proof is the explicit formula of the collectors, which only holds in 
dimension one under a strong diffusivity condition. 

In Ifl3l . an alternative approach is suggested, which is not restricted to the dimension one. 
As in the locally periodic setting, the method is based on a local analysis of the microscopic 
behavior (corresponding to the variable x/e) of the process X s to construct the so-called cor- 
rectors and to identify the limiting process. However, unlike the locally periodic case, these 
correctors turn out to have bad asymptotic properties at a macroscopic scale, in the sense that 
the classical ergodic theory cannot describe their asymptotic behavior. Overcoming this issue is 
the main contribution of fl3l . The main assumption is the uniform ellipticity of the matrix a, 
namely that there exits a constant M > such that for all x,y,X € M. d , 

-^\X\ 2 < (a(u,x,y)X,X) < M\X\ 2 . 

This condition is very convenient for two reasons. From the dynamical angle, it ensures the 
local ergodicity of the process X s . From the technical angle, it provides strong estimates of 
the transition densities of the process X s as well as regularity properties of its generator. The 
control of the process X s , in particular its invariant measure and its tightness, is easily derived 
from this assumption. 
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In this present paper, we intend to improve this latter work by removing the uniform ellip- 
ticity assumption. It is replaced by microscopic ergodicity conditions (Assumption 12.51 ). which 
seem not too far from being minimal to apply classical ergodic theory and then pass to the limit 
in ([]]). The class of considered coefficients then includes possibly degenerate matrices a. In 
other words, we can treat diffusion coefficients a that may reduce to along some directions. 
Under suitable assumptions, we will prove that the process X £ converges to the solution X of a 
SDE with deterministic coefficients, whose generator can be rewritten in divergence form as 




where the so-called homogenized coefficients A and H are respectively symmetric positive 
and antisymmetric. It is worth emphasizing that A may degenerate, even under strong non- 
degeneracy assumptions of the initial diffusion coefficient a. We will prove that the limiting 
diffusion is trapped in a fixed subspace of R d and possesses strong diffusivity properties along 
this subspace. 

We should finally point out that there are only a few papers dealing with possibly degenerate 
diffusion coefficients in the whole literature about probabilistic homogenization of SDEs. In the 
periodic setting, recent advances have been made by Hairer and Pardoux in [0. Their approach 
deeply differs from ours. They allow the diffusion to be strongly degenerate in some area of the 
torus, and even to reduce to over an open domain, provided that the diffusion quickly reaches 
a strongly regularizing area (typically, it satisfies a strong Hormander type condition). Our 
approach does not allow locally such strong degeneracies but does not require any regularizing 
area. As a consequence, we can construct examples that are everywhere degenerate. Moreover, 
the technics used in (5 ) rely on the compactness of the torus and cannot be adapted to the random 
setting. 

The structure of the paper is the following. In section 2, we introduce all the notations and 
assumptions. Our results are stated in Section 4 and an example is given in Section 5. The 
construction of the corrector is carried out in Section 6. Section 7 deals with the regularity 
properties of the process X s such as its invariant measure and the Ito formula. Section 8 is 
devoted to establishing the asymptotic properties of the process X s . Section 9 explains the 
proofs of the homogenization procedure. The tightness of the process X s is treated separately 
in Section 10. 

2 Setup and Assumptions 

Random medium. From now on, d > 1 is a fixed integer. Following Q, we introduce the 
following 
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Definition 2.1. Let (Q,Q,/j,) be a probability space and {t x ;x G M d } a group of measure 
preserving transformations acting ergodically on Q: 

1) VAeg,Vx£ R d , h{t x A) = fi(A), 

2) If for any i£l d r x A = A, then 11(A) = or 1, 

3) For any measurable function g on (Q, Q, p), the function (x, uS) 1— > g[r x u) is measurable 
on {R d x fl, B(R d ) ® Q). 

The expectation with respect to the random medium is denoted by M. Denote by L 2 (Q) the 
space of square integrable functions, by | _ 1 2 the corresponding norm and by (., .)2 the associated 
inner product. The operators denned on L 2 (Q) by T x f(uj) = f(r x u) form a strongly contin- 
uous group of unitary maps in L 2 (Q). For every function / G L 2 (£l), let f(ui,x) = f(r x uj). 
Each function / in L 2 (Q) defines in this way a stationary ergodic random field on M. d . In 
what follows we will use the bold type to denote an element / G L 2 (Q) and the normal type 
f(u, x) (or even f(x)) to distinguish from the associated stationary field. The group possesses 
d generators (throughout this paper, e$ stands for the i-th vector of the canonical basis of Mr) 

(4) Dig = lim he ^ — — if exists, 

h->o h 

which are closed and densely defined. Setting 

(5) C = Span (g*ip;g G L°°(f2), ip G C^°(IR <i )} , with g * (f(u>) = / g{r x uj)ip(x) dx, 

L J J«. d 

the space C is dense in L 2 (J7) and C C Dom(Dj) for all 1 < i < d, with Di(g * ip) = 
—g-kd(p/dxi. Ifg G Dom(Dj), we also have Di(g-kip) = Dig*ip. For / G P|^ =1 Dom(L>j), we 
define the divergence operator Div by Div/ = Yli=i Dif. We distinguish this latter operator 
from the usual divergence operator on M. d denoted by the small type div. 

Locally stationary random fields. Following the notations introduced just above, for a mea- 
surable function / : x R d -> R n , (n > 1), we can consider the associated locally stationary 
random field (x, y) \-> f(r x u, y) = f(oj, x, y) (or even f(x, y)). 

Structure of the coefficients. The coefficients a : x M d -> M. dxd , H : x R d -> R dxd , 
(j \ — ► M dx ' i and 1/ : M. d — > K denote measurable functions with respect to the underlying 
product (j-fields. As explained above, a and H define locally stationary random fields and a 
a stationary random field. H is antisymmetric. We define two new matrix-valued functions by 
a = era* and a = aa*. Furthermore, for some positive constant A, the coefficients a, H, a 
and V satisfy 

Assumption 2.2. (Regularity). For each fixed uj G the coefficients a (u, ., .), H(oj,.,.) 
and a(co, .) are two times continuously differentiable with respect to each variable and are, as 
well as their derivatives up to order two, k-Lipschitzian and bounded by A. V is three times 
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continuously dijferentiable and is, as well as its derivatives up to order three, bounded by A and 
A-Lipschitzian. 

Let us now describe the degeneracies of the matrix a. Roughly speaking, the degeneracies 
of a are assumed to be controlled by the reference matrix a. To be more explicit, let us first 
introduce the 

Definition 2.3. Given a dx d matrix-valued function g :M. d — > M dxa! , a dx d symmetric matrix 
A and a real C > 0, g is said to be (C, A)-controlled ifVy, y' G M. d 

\g(y)\ < CA, and \g(y) - g(y')\ < CA\y - y'\, 

where \M\ = (MM*) 1 / 2 stands for the absolute value of the matrix M (given 2 symmetric 
matrices A, B, the relation A < B means that the matrix B — A is symmetric positive). 

We now precise the control of a by a: 

Assumption 2.4. (Control). We assume that 

M~ l a(uj) < a(u,y) < Ma(uj) 

for some strictly positive constant M and for every (uj, y) G 17 x IR d . Moreover, for any i, j G 
{1, . . . , d} and (u>, y) G x Mr, the matrices d yi a(u>, y), dy iy .a(u>, y), H(uj, y), d yi H(u>, y), 
dy. y .H(uj, y) are (M, a(uj))-controlled. We further assume that 

\a{oj,y + h) - a(uj,y))\ 2 < Ma{u))\h\ 2 

for any y,h G M d and that J^ d e~ 2V ^ dy = 1. 

To ensure the local ergodicity of the process X s , we make the following assumption: 

Assumption 2.5 (Ergodicity). Let us consider the Friedrich extension ( see Mlp- 53] or Section 
5) of the symmetric operator S defined on C C L 2 (Q) by S = (1/2) 2~2ij=x Di(dijDj). This 
extension, still denoted S, is self-adjoint. We then assume that the semi-group generated by S 
is ergodic, that is its invariant functions are \x almost surely constant (see e.g. Rhodes M2\l ). 

Remark. Assumptions \2.2\ may appear restrictive and can surely be relaxed (see 43/ for results in this 
direction in the context of quasilinear PDEs). In particular, the statement of the homogenization property 
only involves the derivatives of order one with respect to y £ K d (see Theorem \3.1h However, it avoids 
dealing with heavy regularizing procedures that are not the purpose of this work. 

Diffusion in a locally ergodic environment. For j = 1, . . . , d, we define the coefficients 

(6) 

I d e 2V%) d 

b j(u,y) = -J^A(a + H)ij(uj,y), Cj(u,y) = -^—J2 d y t { e ~ 2V [ a + H ]ij)( UJ ^y)- 

i=l i=l 
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From Assumption |Z2j the functions bj(ui, ., .) and Cj(u>, ., .) are Lipschitzian so that, for a stall- 
ing point x G K d and e > 0, we can consider the strong solution X s of the following Stochastic 
Differential Equation (SDE) with locally stationary coefficients: 

(7) Xl = x + - f b (X £ r , Xf.) dr + [ c (X £ r , X £ r ) dr + f a (~X £ r , X £ r ) dB r , 

£ Jo Jo Jo 

where we have set X\ = Xf je and B is a standard d-dimensional Brownian motion (the random 
medium and the Brownian motion are independent). We point out that the generator of this 
diffusion could be formally written in divergence form as 

(8) C £ = \e^ £ ^(e-^a + H](u, X /e, X )^). 

i,j=l 1 3 

Notations. For the sake of simplicity, we indicate the starting point x of X s by writing, when neces- 
sary, P| (and W x for the corresponding expectation), this avoids heavy notations as X £ ' x . We can then 
consider the probability measure P £ = M J Rd F E x [.]e^ 2V ^ dx and its expectation W . 



3 Main Results 

Let us now state the main result of this paper. Under the previous assumptions, we can prove 

Theorem 3.1. Homogenization. The law P e of the process X s weakly converges in C( [0, T] ; M. d ) 
towards the law of the process X that solves the following SDE with deterministic coefficients 
( they do not depend on the medium 0, ): 

(9) X t = x+ [ B(X r )dr+ [ A 1/2 (X r )dB r . 

Jo Jo 

The coefficients A and B are of class C 2 and are defined, for y G M. d , by 
(10a) A{y) = lim M[(J + Du x )*a(I + Du x )(.,y)}, 

A — *-0 

(10b) H{y) = lim M[(7 + Du\)*H(I + Du\)(., y)], 

X — ^0 

(10c) B(y) = (l/2)e 2V ^d y (e~ 2V [A + H])(y). 

Formally speaking, for each y £ M. d and A > 0, the entries (u^(.,y)) 1< . < , of the function 
U\(.,y) : fl — ► M. d solve the following so-called auxiliary problems, which are stated on the 
random medium 

V) ~ 2 Yl D i + H jk)Dku\{.,y)] = bi(.,y). 
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Remark. A rigorous description of U\(.,y) is given in Section® In particular, in this degenerate 
framework, the "gradients" Du\ do not exist but along the direction a, that is the only expression 
<x* Du\ can be given a rigorous sense. Because of the control of a and H by a (Assmnption \2.4\) , it then 
makes sense to consider formulae ( llQab and ( II Obi ) (see Section\6\for further details). 

Since the diffusion coefficient a is allowed to degenerate, the reader may wonder whether 
the homogenized diffusion coefficient may also degenerate. The following proposition details 
the structure of the limiting diffusion coefficient A: 

Proposition 3.2. Geometry of the homogenized coefficients. The kernel K = Ker(^4(y)) of 
A(y) does not depend on the point y G M. d where it is computed. For each y G M. d , B(y) G K 1 - 
(K 1 - is the orthogonal complement to K) and there exists a constant QfJ^] ^ ®> sucn tnat 



In other words, for each starting point x G M. d , the limiting process X ( see © ) can be seen as 
the solution of a SDE defined on x + K 1 - with a uniformly elliptic diffusion matrix A. 

4 Example 

Let us consider a simple example in the two dimensional 2-7r-periodic case. The 2-dimensional 
torus T 2 is seen as the random medium equipped with the induced Lebesgue measure, still 
denoted by p, to stick with the notations of the paper . We aim at constructing a degenerate 
homogenized coefficient. For this purpose, let us first define 



where c 7r<Q is a constant, and a = a a* . Choose now any smooth function U : K 2 x M 2 —> 
]R 2x2 , with bounded derivatives up to order 2, 27r-periodic with respect to its first argument 



V(x,y) £l 2 x M 2 , ilrtd < UU*(x,y) < Mid. 

Define V(x,y) G M 2 x E 2 , V(y) = e"^ 2 /vr, a(x,y) = aU{x,y) and H = 0. Let us check 
that these coefficients satisfy all our assumptions. From the smoothness of the coefficients, it is 
plain to see that Assumptions 12.41 and 12.21 are fulfilled. Assumption [23] results from the Weyl 
equipartition theorem (c -zKQ). Theorem B.ll thus holds. 

Let us now prove that A is degenerate and does not trivially reduce to 0. Let us denote by A 
the homogenized coefficient associated to a. From the proof of Proposition I3.21 for any y G M? 
and X G M 2 , we have 




x\ 2 (x,A(y)x) < ogj$x 




x G M 2 and satisfying 



C-\X, AX) < (X,A(y)X) < C(X,AX) = 0. 
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So we just have to compute A. Since <x is constant, it is straightforward to check that A actually 
matches aa* with the help of (|45T ). Indeed, for a given smooth function (p denned on T 2 and 
x £ M 2 , the right-hand side of (T45T ) expands as 

M[\a(Dcp + x)\ 2 } =M[\a*Dip\ 2 } + 2 (5*x, ?*M[D<^]) + (a*x, a*x) 
=M[\a*D i p\ 2 ] + {a*x,a*x) . 

The infimum is then clearly reached for ip = 0. 

Finally, we let the reader check that A = aa* does not reduce to and that the vector 
X K = [1 - c]* satisfies = 0. 

In a general way, because of the various geometries of random media, it is not clear whether 
A is degenerate or not. The reader may find in |3] examples (in a slightly different framework) 
where the diffusion matrix reduces to though the diffusion coefficient a is elliptic over a set 
of full Lebesgue measure, and conversely, an example where <x degenerates and A is uniformly 
elliptic. □ 



5 Construction of unbounded operators 

Throughout this paper, we will need to construct suitable extensions of unbounded operators 
defined on a dense subspace of a given L 2 -space. This construction is always the same and 
follows (31 Ch. 3, Sect 3.] or J9j Ch. 1, Sect 2.], to which the reader is referred for further details 
than those given below. That is the reason why we explain it in a generic way. We also point out 
that the Friedrich extension of S (see Assumption 12.51 ) corresponds to this construction. 

Consider a probability space Q equipped with a probability measure P, a dense subspace 
V of L 2 (Q;F), a positive symmetric bilinear form (•, •) defined on T> x T> (|| • || denotes the 
corresponding semi-norm) and a bilinear form B on V x V that satisfies for any ip,i/j ET> 

(11) a" 1 |M| 2 < Bfatp), B{y^)<a\\y\\U\\ 

for some positive constant a > 0. Let us denote (•, -)2 the canonical inner product on L 2 (Q; P). 

From now on, we will say that the unbounded operator L on L 2 (il;P) is constructed from 
(O, P, (•,•), B) if it is constructed as follows. We consider the inner product II on V x V defined 
by 

and the closure H of V with respect to the corresponding norm. For each A > 0, the bilinear 
form B\ is defined on V x V by 

B x (cp,^) = \(tp^) 2 + B(ip,ip). 

From (TTTb . B\ obviously extends to HI x H (this extension is still denoted by B\). Furthermore, 
it is continuous and coercive on EI x EL Thus it defines a resolvent operator G\ : L 2 (Q,, P) — > EI, 
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which is one-to-one. We can then define L as A — G^ 1 with domain Dom(L) = G\(L 2 (£l, P)). 
This definition does not depend on A > 0. It is readily seen that a function ip G EI belongs to 
Dom(L) if and only if the map ip G H i— ► B\(ip,ip) is L 2 (Q,P) continuous. In this case, we 
can find / G L 2 (f2,P) such that B\(ip, •) = (/, -)2- Then L</? exactly matches / — \ip. Note 
that B(ip, ip) = —(Lip, Tp) 2 for any ip G Dom(L) and We point out that the unbounded 

operator L is closed and densely defined. Moreover, its adjoint operator L* in L 2 (f2;P) coin- 
cides with the operator constructed from (O, P, (•,•), B), where the bilinear form B is defined 
on V x V by B(ip, ip) = B(ip, ip). As a consequence (L*)* = L. 

Notations. In what follows, the notation (FJ, L, Dom(L), (Ga)a>o) = S((S~i,P, {-, •) , i?)) 
means that H, L, Dom(L), (Ga)a>o are constructed from (ft, P, (•, •) , 1?) as explained above. 

6 Auxiliary Problems 

Sefap and notations. Let us now focus on the different operators induced on the random 
medium £1 by the matrices a(-,y) and H(-,y), for each y G M. d . We aim at extending the 
following operators defined on C by 

(12) S^ = - AK-0,y)£j), V> = - D i ((a + H) ij (.,y)D j ), 

i,j=l i,j=l 

according to the method detailed in Section [5] 

The positive symmetric bilinear form (•, - )i * s defined on C x C by 

(13) (^,-0)1 = -(<p,Si/>) 2 = (l/2)(aDcp,Dip) 2 , 

and the associated seminorm || • ||i by ||v?|| 2 = (<p, <p)\. 

For any cp,ip G C, we define the bilinear forms (y is fixed) 

£ 5 (<^) ee -(SW)2 = {ll2){a{',y)Dtp,D1>) % , 
B L (cp, ip) = -{Ly<p, ip) 2 = (l/2)((a + H){;y)Dcp, Dip) 2 . 

From Assumption l2.4l and the antisymmetry of H, it is readily seen that M _1 < B s (ip, (p) 
(resp. M" 1 !!^!! 2 < i3 i (v J ,^))ande 5 (^,-0) < M^y^Hi (resp. B L (<^,^) < 2M|M| 1 ||V||i). 
We can then define 

(Hi, S y , Dom(^), (Gf ) A >o) =E(n, fM, (; -)i,B S ), 
(Hi.i", Dom(^), (Gf )a>o) =S(n, //, (•, Oi.^)- 

Let us additionally denote by (L y )* the adjoint operator of L y in L 2 (f2). Note that S y is self- 
adjoint. 
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We define the space D as the closure in (L 2 (Q)) d of the set {a* Dip; ip G C}. We point 
out that, whenever p,tp belong to C, 2(p,tp)i = (a* Dtp, a* Dtp) 2, so that the application 
G : C — > D, p ^ a* Dp can be extended to the whole space Hi. For each function / G Hi, 
we will note V CT / for Q(f) and this represents in a way the gradient of the function / along 
the direction <r. Similarly, for each fixed y G M. d , we define for any p G Hi the gradient along 
the direction <r(-, y). It will be denoted by V a ^ ,y ^p and is equal to er(-, y)*Dp for any p G C. 
From Assumption I2.4L for each p G Hi, the mapping i/Gl^w SJ al ^ y )p g D is continuous: 

(14) V(y,/i) G (M d ) 2 , \V ai -' y+h) p-V a{ -' y) p>\ 2 2 < M\h\ 2 \\p\\l 
For y G M d and p, ip G C, we derive from Assumption 12.41 

(15) (IV,^2 = -^( J Dv,(a + Ja r )(-,y) J DV) 2 < C|V<V| 2 |V<^| 2 , 

so that we can define a bilinear form T y on the whole space DxD such that \/p, ijj G C 

(16) - (L y p, VO2 = T y (X7 a p, V>). 

Thanks to Assumption 12.21 we can consider the differential dT y of T y defined, for ip, tp G C, 
by dT y (<p,ib) = d y (T y (ip,ip)). From Assumption f2A\ and similarly to ([HI), dT y extends to 
DxD. From Assumption E31 it is then plain to see that the relation d y (T y ($, £)) = dT y (£, C) 
still holds for ^,(£D. 

Whenever a function b satisfies the property: 

(17) Booypec, (&,¥>) 2 <c|M|i, 

we will say that b G H_i and we will define ||b||_i as the smallest constant C satisfying this 
property. 

Solvability and regularity of the resolvent equation. For h G L 2 (VL), u\(.,y) = G^ y h 
belongs to Hi n Dom(L y ) and satisfies \u\(-,y) — L y u\{-,y) = h. Suppose that the right- 
hand side h = h(-, y) depends on the parameter y G M d . We now investigate the y-regularity of 
u x (-, y) from the regularity of y ^ h(-, y) with respect to the norms | • | 2 and || • ||_i. We claim 

Proposition 6.1. Let us consider h : y G R d i-> h(., y) G L 2 (Q) and f : // ; ■ ■ /(., y) G 
L 2 (J7) n H_i. Suppose that there exist C 2 , C_i such that: 

1) the application y 1— > h(.,y) G L 2 (J7) is two times continuously differ entiable in L 2 (Q). 
The derivatives up to order 2 are bounded by C 2 in L 2 (ft) and are C2-Lipschitz in L 2 (Tt). 

2) the application y 1— ► /(., y) G L 2 {Q) n H_i is two times continuously differentiable in 
H_i. The derivatives up to order 2 are bounded by C_i in H_i and are C -\-Lipschitz in H_i. 

Then, for any A > 0, the solution u\(.,y) G Hi n Dom(L y ) of the equation 

(18) Xu x {., y) - L y u x (.,y) = h(., y) + /(., y) 
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is two times continuously differentiable in Hi with respect to the parameter y G M. d . Further- 
more there exists a constant D^J\ > 0, which only depends on M, C_i, such that the functions 
9x(-,y) = u\(-,y), dyux{.,yld% y u\(.,y) satisfy the property: V(y,h) G R 2 , 

(19a) M9x(,y)\l + \\9x(;y)\\i < + C 2 2 A), 

(19b) X\g x (.,y + h)- g x (.,y)\ 2 2 + \\g x (.,y + h) - g x (.,y)\\l < Lfcjfl + C%/\)\h\ 2 . 

Proof: The proof is readily adapted from [ 13 , Prop. 4. 1]. The method consists in differentiating 
the resolvent equation (fT8l) with respect to the parameter y G M. d . In the uniformly elliptic setup 
|[T3l Prop. 4.1], this can be earned out thanks to the differentiability and the boundedness of 
a, H and their derivatives up to order 2. In the degenerate setup, we need to control the matrices 
a and H, as well as their derivatives up to order 2 with respect to the parameter y G M. d , by the 
matrix a (see Assumption 12.41 ) in order to differentiate the function y t— > u x (-, y) in Hi. □ 

Auxiliary problems: construction of the correctors. The end of this section is now devoted 
to the study of the solutions of the so-called auxiliary problems, that means the solutions u\(.,y) 
(i = 1, . . . , d) of the resolvent equations 

(20) Xu\(.,y)-Lyui(.,y) = b i (.,y), 

where bi(.,y) = (1/2) Ylj=i Dj [( a + H)ji{-> y)] ■ The weak form of the resolvent equation 
then reads for (p G C 

(21) X(u\(.,y),<p) 2 + Ty(^u\(.,y),V^if) = -(l/2){(a + H)(.,y)e i ,D V >) 2 . 
Having in mind to apply Proposition 16. 1[ we first prove 

Lemma 6.2. The mapping y bj(., y) G L 2 (f2) PlH_i is two times continuously differentiable 
in H_i, and the derivatives are bounded and Lipschitzian in Hi. 

Proof: First note that for each cp G C, 

(bi(.,y),<p) 2 = -(l/2)((a + H)(.,y)e i ,D l p) r 

From Assumption 12.41 we easily deduce that bi(.,y) G HL_i and that the mapping y G M d i— > 
bi(.,y) G H_i is bounded and Lipschitzian. 

From Assumption 12.41 again, it is readily seen that the H_i derivatives of foj coincide, for 
1 < k < d, with the classical derivatives d Vk hi and 

(d yk bi(;y)^h = -(l/2){(d yk a + d yk H)(.,y)e l ,D ( p) 2 < CMl 

Since d yk a(uj) and d yk H(uj) are (M, a (lo)) -controlled, the derivatives are bounded and Lips- 
chitzian in Hi. The same job can be carried out for the second order derivatives. Details are left 
to the reader. □ 
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From Proposition 16. II (with h = and / = b{), the mapping y i— ► u\(.,y) is two times 
continuously differentiable in Hi. We now investigate the asymptotic behavior of u\ as well as 
its derivatives, as A goes to zero. 

Proposition 6.3. For each fixed y G M. d and 1 < i < d, the family (V CT, u A (., converges to a 
limit £i(.,y) G L 2 {Q) d as A goes to 0. 77ie s« property holds for the derivatives, namely that 
the families (V a d Vj u\)\, (V 7 dy yk u\)\ (1 < i,j, k < d) respectively converge to <9y 3 £j(-, y), 

dy jy . k $,i{-,y) in L 2 (Q) d . Furthermore, we have 

AK(., y)\ 2 2 + X\d yj u{(.,y)\l + \\d$ m u\(.,y)\2 0, as A feiufc to 0, 
a*zc?, each function g A (., y) = u\(.,y),d yh u\(. } y),d ykyi u x (., y) satisfies the property: 

(22) M9x(,y)\l + \\9x(;y)\\l<q5j] 

(23) A| ffA (., !, + /»)- g x (., y)\l + \\g x {., y + h) - g x (.,y)\\l < 

for every y,h G R d , where Q^j\ is a positive constant independent of A > <3«<f y G M rf . 

Proof: The proof does not deeply differ from Proposition 4.3 in lfT3"l . but we nevertheless 
set it out because of its importance. From dl9a| ) (note that C 2 = 0), we get A|it A (.,y)|2 + 
\V°u\(., 2/) || < C. Denote by y) G L 2 (ft) d a weak limit of the family (V ? w\(., y)) A as 
A goes to 0. Passing to the limit in (|2~TI ). it is plain to see that \/ip G C 

(24) T«'(e i (.,y),V SF v) = -(l/2)((o + lT)(.,y)e i ,D ¥ >) 2 . 

Since is coercive on D x B, this proves the uniqueness of the weak limit in B. Gathering 
d2Tb and d24]), we get 

(25) X(u{(.,y), <p) 2 + ^(VXf^j.VV) = T»&(., y), V*<p). 
Choosing u x (.,y) = cp yields: 

\\u\(., y )\ 2 2 +Ty(V°u\(.,y),V°u\(., y )) < T"fe(.,y),^(.,y)) + e(A), 

where the function e(A) exactly matches T^f^., y), V CT w A (., y) — y)) and thus converges 
to as A goes to 0. Hence limsup A ^ T^(V ? u A (., y), V 5 w A (., y)) < T» (£,(., y),^(., y)). 
Denote by T 5 the symmetric part of 

T 5 (vp, ^) = (1/2) [T y (<^, ^) + T y (</>, cp)) , ^£D. 

From Assumption 12.41 and the antisymmetry of H, we have 

Ar l {a*D(p,a*D<p) 2 < T s (V°cp, V B (p) < M(a*Dcp,a*D<p) 2 , cpeC. 
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By density arguments, the quadratic form associated to T defines a norm on D equivalent to the 
canonical inner product. Moreover, we have just proved that the family (V CT tt\(., y))\ is weakly 
convergentinDtol^.^jandlimsup^oT 5 ^^^,^^^^,^) < T s fy., y),^, y)) 
Thus the convergence is strong with respect to the norm on D associated to T , and conse- 
quently (V CT, u^(., y))\ strongly converges in (L 2 (Q.)) d to ^(., y). From this together with (1251) . 
we get 

A|ul(.,y)H + | W A (.,y) - - as A - 0. 

This proves the first part of the statement for the function u\(.,y). The second part results from 
Proposition 16. 1[ statements (1 19ab and (1 19bb (with C2 = 0). The same job can be carried out for 
the successive derivatives of u\(., y) up to order 2. □ 



7 Dynamics of the process X £ . Preliminary results 

Notations. All the results of this section are valid for any value of the parameter e. However, to 
simplify the notations, we choose e = 1 and thus remove the parameter efrom the notations. So 
the process X stands for the process X £ defined by (Q. Finally we denote by Fy the probability 
measure e~ 2V ^ dy ® dfi on Q x M. d and by My the coresponding expectation. 

This section is devoted to the study of the SI x Revalued process (txuj,X), such as its 
invariant distribution and the Ito formula. Since these properties are more easily established 
when the process X possesses regularizing properties, namely that the diffusion coefficient a 
is uniformly elliptic, most of the following proofs are carried out through vanishing viscosity 
methods, that is, in considering a family of non-degenerate diffusion processes that converges 
to X. 

Invariant distribution. Let us introduce a standard d-dimensional Brownian motion B in- 
dependent of B. For each fixed (oj, n) G Q X M* and for any x G M. d , we define the Ito process 
X n as the solution of the SDE (with the convention n" 1 = if n = 00) 

X? = x + [\b + c - n- l d y V){u, X?, X") dr + f a(to, X?,X?) dB r + (n/2)- 1 / 2 ^. 
Jo Jo 

Note that, for n = 00, X°° coincides with the process X. For n £ IN*, the process X n 
defines a continuous semigroup P n on C{,(M d ) (continuous bounded functions). Its generator 
C n coincides on C 2 (R d ) with 

(26) C n = l -e 2V ^ Y, d ^ {e- 2V[x \a + H + n^ld)^, x, x)d Xj ■ ) . 

hi 

For n G IN*, it is well-known that the distribution of X™ (t > 0) admits a density p n (u!,t,x,-) 
with respect to the Lebesgue measure (cf. [14, Sect. II.2]), which is bounded from above 
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by a constant C that only depends on A, n, t. Thus the semigroup associated to X n (n G 
W) continuously extends to L 2 (R d , e~ 2V ^ dx). Let us denote by (£ n )* the adjoint of C n in 
L 2 (M. d , e~ 2V ^ dx), which coincides on C 2 (R d ) with 

(27) (£")* = ^e 2V ^^d Xi (e- 2V ^ x Ha- H + n- 1 Id) ij {u;,x,x)d Xj •)■ 

1,3 

Now, for (p,ifi G C c °°(IR d ), let us compute / Rd £ n P t n v?(x)V'(x)e- 2y ^ dx. From d, 
-P/V G C 2 (R d ) so that C n Ppip can be computed with the help of ( |26l ). By integrating by 
parts, we obtain 

(28) / £ n P?<p{x)^{x)e- 2V ^ dx = [ P^{x){C n )*^{x)e~ 2V{ - x Ux. 

Moreover, we have C n Pp<p = P/'X'V G C 6 (R d ). Choose now a function £» G C c °°(]R d ) that 
matches 1 over the ball B(0; 1). Define i^rn(x) = g(x/m). It is readily seen that the sequence 
(£ n ip m ) m is bounded in L°°(]R rf ) and uniformly converges to on the compact subsets of R d . 
Thus, choosing ip = ip m in (|28T ). and passing to the limit as m goes to oo, we get 

(29) V^GC c °°(M d ), / £ n P t n (^(x)e- 2y(a;) dx = 0. 

In particular, for any cp G C^°(R d ), f Rd P[ l (p(x)e~ 2V ^ dx = f Rd Lp(x)e~ 2V ^ dx, in such a 
way that, by density arguments, the probability measure e~ 2V ^ dx is invariant for the process 
X n (n > 1). Then classical arguments of SDE theory ensure that the sequence of processes 
{X n ) n converges in law in C([0,T];M. d ) to the process X as n goes to oo. We deduce that 
J Rd Pt(p(x)e~ 2V ( x ^ dx = J Rd ip(x)e~ 2V ^ dx holds for ip G Cb(M. d ). The semigroup associated 
to X thus extends to L p (M d ; e~ 2V ^ dx) for p > 1 and the probability measure e~ 2V ^ dx is 
also invariant for this semigroup. 

Finally, for each ip G Cb(Q x M. d ) (i.e. for each fixed u G f2, the function x i— > (^(r^u;, x) 
is continuous and bounded by a constant independent of lo) and n > 0, we deduce from the 
previous remarks and the invariance of the measure fi under space translations that 

(30) E[(p(Txnu,X?)] = M v [cp(r x u;,x)} = Myfo>(w,x)], 

so that the mapping (p G Cj,(f2 x M. d ) i— > Pp(ip) = E, x [ip(Tx n uJ, X™)\ continuously extends to 
LP(0 x R d ; F v ) for any p > 1 and OU]) holds for ip G X R d ; F v ). 

Ito's formula. We now aim at establishing the Ito formula to the process (txui,X) and 
to the function (x,y) h-> u\(ui,x,y), where u\ is the solution of the resolvent equation (PT8T ). 
with functions h(.,y) and f(.,y) satisfying the assumptions of Proposition I67TI This latter 
proposition describes the regularity of u\ with respect to the variable y. Due to the possible 
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degeneracies of a, the difficulty actually lies in the regularity with respect to the parameter 
x G R. d . To apply the Ito formula and get round technical difficulties, we use viscosity methods 
again, namely that we look at the operator A — L y — re -1 A for n G IN*. Obviously, there is no 
difficulty in solving the corresponding resolvent equation with the techniques used in Section [6] 
(it suffices to replace a by a + ?7, _1 Id and to choose a = Id) 

(31) A«! n) (.,y) - (p + n- 1 A)4 n) (-,y) = h(-,y) + f(.,y). 

The strategy then consists in applying the Ito formula in the non-zero viscosity setting and 
then in letting n tend to oo. Thanks to the regularizing parameter n G IN* , the Ito formula holds 
in the non-zero viscosity setting (cf [13 , Sect. 5]). The following formula thus holds 

du^(X?,X?) =(Auf -h-f){X?,X?)dt + [c-n- 1 d y V]-Du i Z\x?,X?)dt 

(32) + (V^it^VC*?, X t n ) dB t + n-WiDvlpHX?, X?) dB t 
+ bdyU^iX^X?) dt+[c- n^dyV] ■ d y uf\x?,X?) dt 
+ {d y u^Ta(X?,X?)dB t + n- 1 l\d y u ( ?)(X?,X?)dB t 

+ (l/2)trace([a + n^Idjd^u^X^ , X t n ) dt 
+ trace([a + n^I^Ddyuf*)^, X?) dt. 

Having in mind to let n tend to oo in (l32l ). let us now describe the behavior of w" as n tends to 
oo. We first claim: 



Proposition 7.1. 

(33) lim \\u^\.,y)-u x (.,y)\ 2 + \\u^\.,y)-u x (.,y)\\ 1 +n- 1 \DuP(.,y)\l 

n— >oo L 

and that there exists a constant D^^( independent ofn and y G M. d ) such that 

(34) \u[ n \;y + ft) - + Huj^.y + fc) - «?>(■, v)||? 



0. 



+ n- 1 |D«^(-, 2/ + ft) - < Z^Jftl 2 . 



Moreover, the same properties hold for the sequences {d yk u x ) n , (dy kyi u \)n an d their cor- 
responding limits (d yk u\) n , {dy kyi u\) n , for 1 < k, I < d. 

Proof. Since the proofs of d33l and d34l ) can be adapted from the proof of Proposition |6.3[ we 
just set out the guiding line of (l33l) . 

To clarify the notations, we forget for a while the dependence on the parameter y. First 
multiply (|3TT) by and integrate with respect to the measure fi so as to obtain the estimate: 

X\uP\ 2 2 + |V 5 4 n) |2 + n'^Duf^l < C 
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for some constant C only depending on |h||/A and From this estimate, we deduce 

that the family (n~ l Du^) n strongly converges to in (L 2 (Q)) d as n — > oo and that, up to 

extracting a subsequence, the family (u^) n weakly converges in Hi as n — > oo. Multiply once 
again d3TT ) by a test function ip G C, integrate with respect to the measure \i and then pass to the 
limit as n — > oo to identity the weak limit in Hi as being necessarily equal to u x . So the whole 
family (u^) n is weakly convergent in Hi (not up to a subsequence). It just remains to prove 
that the convergence actually holds in the strong sense. We can integrate ( 1311 ) and ( fT8l ) against a 
test function ip G C. Since the right-hand sides of (I3TT) and (fT8l) coincide, this yields: 

A(4 n) , <p) 2 + T«(V'u^,V s <p) + n- l (Du^,Dcp) 2 = \(u x , <p) 2 + T^V^ Al V>). 

(n) 

Choose cp = u x and pass to the limit as n — > oo and get 
Jim (A| W f|i + T^(V^ ) , V 5 4 n) ) + n- 1 |Z?^ ) ||) = A|u A |g + T«'(V af u A> V'u A ). 

As in Proposition 16.31 this is sufficient to establish the strong convergence of (u^) n in Hi and, 

consequently, the convergence n~ 1 \Du^\2 — > as n — > oo. □ 
We are now in position to conclude. Going through formula (l32l) . we are faced with func- 
tional of type g n (X™ , X™) dr (concerning the martingale terms, it suffices to work on their 
quadratic variations), where My [\g n — g \] — > as n tends to oo and 

(35) V(y, h) G R d x M d , \g n (.,y + h) - g n (.,y)\ 2 < C\h\ 

where the constant C depends neither on n G IN" nor y, h G M. d . From Lemma 1731 below, 
we prove the convergence of the functional towards go(X r ,X r ) dr in IP-probability and as a 
consequence the 

Theorem 7.2. Let h, f be two functions satisfying the assumptions of Proposition \6J\ Let u x 
be the solution of the resolvent equation: 

\u x (;y) - Ly Ux (-,y) = h(-,y) + f(-,y). 

Then the following ltd formula holds (we reintroduce the parameter s): 

edu x (X £ t , XI) =e~ 1 (\u x -h- f){T t ,X E t ) dt + c- Du x (T t , X £ t ) dt 

+ (v^ Ux y(x £ t , x!) dB t + bd yUx (x e t , x!) dt 

+ e(d y u x )*a(X £ t ,X e ) dB t + ec ■ d y u x (X £ , X £ t ) dt 

+ (e/2)trace(a# u x )(X £ , X s ) dt + tT&ce(aDd y u x )(X e , X £ t ) dt. 
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Lemma 7.3. Consider a sequence of functions g n G L'(!l x M rf ;Py) (n > 0) smc/i 
My[|f/ n -0 o |] ^ Q as n ^ oo and for any (y,h) G M d x M d , |g n (., y + /i) - gr n (., y)| 2 < C|/t| 
/or some constant C that depends neither on n nor y,h G 
TTwrn E[| 9n (X™, X?) - g (X r ,X r )\] -»0asn->0. 

Proof: First, suppose that g is bounded. Let us consider a smooth mollifier p : R rf — > K and 
o G C^°(R d ) such that o = 1 over the ball B(0; 1). We define for m, q > 1, p m (-) = m d p{m •), 
Q q {-) = e(-/o) and 9™' 9 (u, x ) = Jr<* 9o{t-x'U, x')g q (x')p m (x - x 1 ) dx' . Then, from OOJl, 

-go(X r ,X r )\] <M v [\g n -g \]+2M v [\g^ q -g \] 

+ E[\g™' q (X?,X?)-g™' q (X r ,X r )\}. 

With classical convolution techniques, we can prove that m, q can be chosen large enough to 
make the term 2My — gr |] small. Then, from the Lipschitz regularity of the coefficients 
(Assumption 12.2b . the classical theory of SDEs ensures that E x [sup 0<i<T \X£ — X t \ 2 ] < n~ 1 D 
for some constant D that only depends on M, A and T. For each fixed m, q > 1 and wed, the 
function x \— ► g™' q (x, x) is continuous with compact support so that J Rd E x [\g™' q (X™, X") — 
g™' q (X r , X r )\]e~ 2V ^ dx — > as n — > oo. Then, the Lebesgue theorem ( y™' 9 is bounded 
independently from cj ) proves that K[\g™' q (X™, X™) — g™ ,q (X r , X r )\] converges to as n goes 
to oo. Therefore, n can be chosen large enough to make this latter term small. Finally, from 
the assumptions of the lemma, even if it means considering larger n, the term My[|g n — g \] is 
small too. The proof is then easily completed in the case when g is bounded. 

If g is not bounded, it suffices to consider for n > and R > 0, g% = max(— R; min(<jf n ; R)). 
It is readily checked that the sequence {g^)n still satisfies all the assumptions of the lemma in 
such a way that 1L[\g*(X?, X") - g(?(X r ,X r )\] -> as n 0, for each fixed R > 0. Then, 
from m, E[\g*(X?,X?) - g n (X?,X?)\] < M v [\g% - g n \] and 

Km lim M v [\g% - g n \] = lim M v [\g* - g \] = 0. 

R->oo n— >oo R—ioo 

Since we have 

n\g n (X?,X?)-g (X?,X?)\}<E[\gX(X?,X?)-g n (X?,X?)\] 

+ E[\g*(X?,X?) - gK(X r ,X r )\}+Ei\gX(X r ,X r ) - go(X?,X?)\], 

the proof is then easily completed in this case too. □ 

8 Asymptotic Theorems 

Classical ergodic theorem. In this section, we aim at exploiting the asymptotic properties of 
the process X s , more precisely Assumption 12.51 in order to describe the asymptotic behavior 
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of functionals of type L ^>(X £ r ,X £ .) dr for a suitable locally stationary random field >f. The 
classical ergodic theory leads us to guess that the local ergodicity assumption 12.51 makes this 
functional average with respect to its first variable. More precisely, 

Theorem 8.1. (Ergodic Theorem) Let us consider * : OxI^ R such that My [\^\] < +00. 
Denoting ^(y) = M.[$f(-,y)], the following convergence holds: 

(36) 1 £ [ sup I / ^(X £ r ,X £ )dr- \ ^{X £ )dr\ 2 } ► 0. 

0<s<t Jo JO e ~*° 

Proof: This result can be proved in the same way as lfT3l Th. 6.1]. The only difference consists 
in establishing: g £ Dom(L y ) C Hi and L y g = implies that g is constant \i almost surely. 
In the uniformly elliptic setting, it turns out that the derivatives Dig reduce to and, as a conse- 
quence, g is constant. In the degenerate framework, we need to use Assumption I2.5l as follows. 
From Assumption E31 ||flf||i < MHffH^ = -(g,L v g) 2 = 0. In particular, B sv (g,-) = 0. 
Hence g £ Dom(S v ) and S y g = 0. Thus g is constant (Assumption I2.5I ). □ 

Asymptotic theorem for highly oscillating functionals. Theorem 18.11 describes the asymp- 
totic behavior of functionals of type J * ^(X 6 ., Xf.) dr in order to pass to the limit in How- 
ever, as explained in lfT3l . additional difficulties arise in the random setting in comparison with 
the periodic one. In particular, we must describe the asymptotic behavior of the functional 
/q ^ £ (X £ r , X £ ) dr for a family (* e ) e that need not be convergent in L l (£l x R d ; Py) but satis- 
fies a sort of uniform Poincare inequality. Unlike |[T3l Theorem 6.3], technical difficulties due 
to the degeneracy of the diffusion coefficient a occur. In particular, because of the lack of Aron- 
son type estimates, the tightness of the process X £ is not obvious. To prove this tightness, all 
asymptotic convergences need be established in C([0, T]; U. d ) (note the sup in (I38T)). This is one 
of the main difficulty of Theorem 18.21 below in comparison with the uniformly elliptic setting 
(see H3l Theorem 6.3]). The strategy consists in expressing f£ ^ £ (X e .,X £ ) dr as the sum of 
two martingales thanks to time reversal arguments, and then in using the Doob inequality. The 
Poincare inequality d3"7T ) ensures that the martingales possess suitable asymptotic properties. 

Theorem 8.2. (Ergodic theorem II) Let us consider, for each e > 0, a function ^f £ E L 2 (fi x 
M d ;Py) satisfying the following Poincare inequality: for any ip(u,x) = x(w)f}(x), (x, f?) € 
C x C™{R d ), 

(37) My[* e¥ >] < C £ (M y [| CT *(7J + edj>| 2 ]) 1/2 , 
for some family (C e ) e >o satisfying eC £ — ► as e — ► 0. Then 

(38) K £ \ sup I f ^ £ (X £ r ,X £ )dr\ 2 } > 0. 

0<s<t Jo 
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Proof: In what follows, we say that ip G Cn if <p(uj,y) = %(u;)^(y), where (%, g) G C x 
C£°(M d ). We aim at constructing, as prescribed in Section |5j the unbounded operators on 
L 2 (J7 x M. d ; Py) that coincide on Cn for n G M* with (here we use the convention n^ 1 = if 

71 = OO) 

(39) S n ' £ ip =(l/2)e 2V ]T (A + e^)[e" 2V '(a + n- 1 Id) lJ p j +^ 3 )^], 

ij=l,...,d 

(40) L n ' s <p =(l/2)e 2V (A + e9 ^) [ e_2F ( a + H + ™~ lld W(^ + • 

ij=l,...,d 

For e > 0, n G M* and (p, i/> G Cn, we define the corresponding bilinear forms 

(41) {cp, ^) rei£ = (1/2)M V [(£>¥> + sd y cp)*{a + n^Id)^ + eS„^)] , 

(42) B rh£ {p,ijj) = (1 /2)M V [(Dtp + sd y <p)* (a + H + n~ l U)(Dt{> + ed y <4>)}. 

Clearly, (•,•)„ £ is positive symmetric (denote by || • || nje the corresponding seminorm). Note 
that, for each fixed e > 0, the seminorms (|| • || n ,e)ne]N* w& a U equivalent. Moreover, for 
n G W, \\p\\l t£ < B n>£ (cp,cp) and B n>£ (cp,iJ>) < 2M 2 ||^|| ni£ ||-0|| ni£ for any cp, ip G C n (see 
Assumption 12.41) . From Section [5j we can define 

(M nte ,S^,Dom(S n > s ),(G S x ' n ' £ ) x>0 ) = H(fi x M d , Py , (•, .) n , e , <•, -) n , e ), 
(M n , £ , L"> £ , Dom(L™' £ ), (G^ n ' e ) A>0 ) = E(fl x R d , Py, (■, , B n , £ ). 

and we denote by (L n ' 6 )* the adjoint operator of L n ' 6 in L 2 (fi x R d ; Py). 

Let us now consider a family (\f e ) £ of functions in L 2 (f2 x Py) satisfying (l37l) for some 
family (C e ) e>0 such that eC e -> as e -> 0. Fix n G M*. Define <£ n>£ = G^ n ' e (* e ), which 
satisfies e 2 M v [cp n ^\ + (c^, -0)^ = M v [& £ ip] for any ^ G H n)E . Choosing i/> = V n>£ , 
using (1371 ) and the standard estimate ab < a 2 / 2 + 6 2 /2 leads to 

e 2 My[|^ i£ | 2 ] + \\<p n Jl £ = M,[* eV „ iE ] < C £ V2\\cp n J 0>£ < V2C £ \\cp n;£ \\ n!£ 

<C? + \\<Pnjle/2 

in such a way that 

(43) e 2 M v [\cp n J 2 ] + \\cp n X,e/2<C 2 . 

Once again, to apply the Ito formula, we use vanishing viscosity methods in order to get 
round the lack of regularity of cp n £ because of the degeneracy of a. In the non-degenerate 
framework (n > 1), from |[T3l Proof of Lemma 6.3], standard convolution technics provide us 
with a H n e -sequence (¥>™ e )meisr of smooth functions, namely that for each fixed uj G $7 the 
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function x i-> <p™ £ (r x/e u}, x) is a C°°(IR d ) -function, such that M v [\tp™ e - <p n J 2 + \S n ' £ (p™ £ - 

S n ' 6( Pn,e\ 2 } + Wfnfi ~ ¥Vlln,e ~> as m § oes to °0- 

We are now going to use a time reversal argument. Let us consider the process (introduced 
in Section Q 

Xt' £ =x+ f {e- l b+c-n- l d y V){uj,X n r ' e ,X?> £ )dr+ f a{uj, IT/, X r r l ' £ ) dB r +{n/2y 1/2 B t , 
Jo Jo 

where X™' £ = X™' 6 /e. As explained in Section |7J its generator coincides on C 2 (R d ) with 

p 2V{x) 

C n,e = _____ £ Q x% ( e -2V(x) {a + H + n^ld)^, x/e, x)d Xj • ) 
i,j 

and admits e~ 2V ^ dx as invariant measure. Furthermore, for a fixed T > 0, the generator of 
the time reversed process t X^_ t with initial law e~ 2V ^ dx coincides with the adjoint of 
C n ' £ in L 2 (R d ; e ~ 2V ^ dx). For each cp G C 2 (R d ), it exactly matches 

{C n ' £ ) V = £ {e- 2V(x \a -H + n-hd^jiu, x/e, x)d Xj <p) 

As a consequence, observe that, for any < s < t < T, 

vl) £ (X ,£ ,Xr) =< £ (Xr,X s n ' e ) + [\c n ' £ (rt] e (-/e,-))](X n /,X?' £ )dr 

J s 

+ (M^ n ' E -M™' n ' £ ), 

where J\4 m ' n ' £ is a martingale with respect to the forward filtration (J r ™' £ )o<t<T and T™' £ is the 
cr-algebra on M. d generated by {X"' £ ; < r < t}. In the same way, 



where _A4 m ' n,e is a martingale with respect to the backward filtration (Gt' £ )o<t<T and Q £ is the 
cr-algebra on R. d generated by {X™' £ ; t < r < T}. Add these two expressions: 

-2e~ 2 f S n > £ <p™ £ (Xr' £ , X^ £ ) dr = (M™' n > £ - M™> n > £ ) + (M™' n ' £ - M™> n > £ ). 

J s 

We further mention that the quadratic variations of both martingales exactly match 
e~ 2 [\(D + ed y )^ e M(D + ed y )^J '(X?' e , X?< £ ) dr, 

J s 
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in such a way that the Doob inequality yields 

E £ \ sup | f S^^ e {X ,e ,Xr)drf\ <WTe 2 \\<p™ £ \\ie- 

L 0<s<i JO J 

Letting m go to oo, reminding that e 2 ip n £ — S n ' £ ip n e = \l/ e and using (03]) leads to 

E £ [ sup | l\ e QT; e ,X^)dr\ 2 ] <32Te 2 \\<p n£ \\l £ + 2Te 4 M v [\<p n r\ 2 }<68Te 2 C 2 . 

0<s<t Jo 

We then complete the proof in letting n go to oo and in using the fact that X n,£ converges in 
C( [0, T] ; R d ) towards X £ as n goes to oo. □ 



9 Proof of Theorem [3J] and Proposition 1372 



Proof of Theorem 13.11 Section [TO] below is devoted to proving the tightness of the family of 
processes (X e ) £ in C([0, T]; W 1 ). It remains to prove that there is a unique possible weak limit 
for all converging subsequences. 

From now on, the corrector u\ (A > and 1 < i < d) stands for the solution of (l20l >. 
Applying the Ito formula (Theorem 17.21 ) to the correctors leads to 

dX £ t = - edu £ 2 (X £ t , XI) + e(d y u £ 2)*a(X £ t , X s ) dB t 

+ [eu £ 2 + ec • d y u £ 2 + (e/2)tTa,ce(ad 2 y u £ 2)](X £ t , Xf) dt 
+ [bd y u £ 2 + c • (J + Du £ 2) + trace(aDd y u £ 2)](X £ , X^) dt 
+ [o + Du £ 2cr}{X £ t ,X £ t )dB t 
=d&l' £ + d@ 2 t ' £ + d@ 3 t > £ + d&* ' £ 

Concerning the first term, we have E £ [|6 t 1,e | 2 ] < (1 + T)e 2 M v [\u £ 2\ 2 + M 2 \d y u e2 \ 2 ] for 
< t < T. This latter quantity converges to as e goes to from Proposition 16.31 The same 
job can be carried out for @ 2,£ and the same conclusion holds. 

The main difficulty actually lies in the term 3 ' e , especially in the part corresponding to 
bd y u e 2. Concerning the remaining part c • (I + Du e i) + tvace(aDd y u e 2), it is readily seen 
(see Proposition [63]) that it converges in L 2 (Q, x M. d ; Fy) and thus Theorem l8.1l can be applied. 
As a consequence, we have 



E £ [ sup | / [c-(I+Du £ 2)+tTace(aDdyU £ 2)}(X £ ,X £ )dr- / <S>(X £ ) dr\ 2 ] -> as e -> 0, 

0<t<T Jo Jo 

where $(y) = lim^^ M[c • (/ + Du\) + tvace(aDd y u\)(-,y)]. It remains to treat the term 



(bd y u e 2) £ . Note that the L 2 -norm of bd y u £ 2 need not be convergent. That is why we have in 
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mind to use Theorem 18.21 Up to introducing new correctors, we will prove that bd y u £ 2 can be 
divided into two parts, satisfying respectively Theorems 18.11 and 18.21 To understand how this 
decomposition occurs, let us consider a test function ip G Cjj. Then two successive integrations 
by parts yield, for 1 < i,j < d, (we use the convention of summation over repeated indices) 

My [bjdy^cp] = (1/2)M V [D p (a + 1 1 ) l , J i), / .u\, i p[ 

= -(l/2)Mv [(a + H) pj {D p dy 3 v! £2 ip + Qy.v^D p (p)] 

= -(1/2)M V [(a + H) pj (D p d yj ui 2 <p + d y ^(D p + ed yp )<p)] 

+(e/2)M v [(a + H) pi d yj u^d Vp ip\ 
= -(1/2)M V [(a + H) pj D p d yj ul 2 <p] - (l/2)My [d y ^{D p + e^)^] 

-(e/2)My [d yp (a + H) pj d y ^<p + (a + JT)^ 2 .^^] 

+eM v [(a + HJp^-^.t*^^ V] . 

So, for 1 < z < d, define the correcting part Corr*(u;,y) = (e/2)d Vp (a + H) p jd yj u l e2 + 
(e/2)(a + H) pj dl.y p u l E2 - e(a + H) pj d yj u l e2 dy p V , the L 2 -converging part Conv|(w,y) = 
-(l/2){a+H) pj D p d yj u l e2 and L 2 -diverging part Div*(w, y) = [&j<9 % .u* 2 +Corr* -Conv*](w, y) 
From the previous calculation, Div* satisfies the "Poincare inequality" (1371) . namely that for 
any function <p in C n , My [¥ e¥ >] < (My[|<9yW* 2 | 2 ]) 1/2 (My[|(Z) + ea y )vj| 2 ]) 1/2 . Moreover, 

Proposition 16.31 ensures that e(My[|(9j,w* 2 | 2 ]) — > as e goes to 0. Consequently, (l38l) holds 
for Div* . Thanks to Proposition [63J the family (Corr* ) e converges in L 2 {n x M d ;Py) towards 
0. As a consequence, W [( J * Corr* (X e r , Xf.) dr) 2 ] tends to as e goes to 0. Then, Theorem 
EH ensures that W [ sup < t < T | /„* Convex 2 , X E r ) dr - jj f(X r £ ) drf] -> as e 0, where 
r(y) = lim^_>o — (l/2)M[(o + H) p jD p d yj u\{- , y)] . To sum up, this proves that 

(44) E e [ sup | I b-d y ui 2 (X £ ,X £ r )dr - [ T(X £ r ) dr\ 2 ] -> 

0<t<T JO JO 

as e tends to 0. 

Concerning the martingale part G 4 ' £ , it suffices to apply Theorem 18. II to the quadratic vari- 
ations. 

Hence each possible limit point X in C[0,t];R") of the process X s must solve the mar- 
tingale problem X t = x + J Q * B(X r ) dr + J * A^ 2 (X r ) dB r , where the entries of B are given 
by 

B~i = lim M[ - (l/2)(a + H) pj D p d yj u{ + c,((% + + a^-D^wl] 

= lim M[(l/2)(o + H) pj D p d yj u{ + + D jU {)] 

A — >\) 
e 2V 

= — d v 3 (e- 2V limM[(a + H)^ + D p u\)] ) 
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Thanks to Proposition 13.21 it is readily seen that the coefficients B and A are two times 
continuously differentiable with bounded derivatives up to order two. In particular, they are 
Lipschitzian and there exists a unique solution to the corresponding martingale problem. □ 

Proof of Proposition 13.21 The strategy consists in introducing the homogenized diffusion co- 
efficient associated to the operator S and in comparing it with A(y). So we define the d x d 
nonnegative symmetric matrix A as the unique symmetric matrix satisfying (this is the classical 
variational formula for the homogenized coefficient associated to S, see iPTOl for further details) 

(45) VxGM d , (x,Ax) = inf M\\a* {Dip + x)\ 2 }. 

Due to Assumption I2.4[ we have for each function ip G C, 

M~ x {x,Ax) < M~ 1 M[\a*(D<p + x)\ 2 ] < M[|tr*(-, y)(D<p + x)\ 2 ] . 

Since C is dense in Hi, we can choose (p = u\(-, y) ■ x and then pass to the limit as A tends to 
0. We obtain M _1 (a;, Ax) < (x,A(y)x). 

Now we turn to the auxiliary problems (subsection [6]). Denoting by L the closure of 
{ct*C, C G L 2 (n- W d )}, we can extend T y to the whole L as follows 

(46) vc, e e L 2 (n, R d ), tv*C °*o) = (1/2) ([a + h](-, y )c, e) r 

The underlying quadratic form is still denoted by T y {-). Furthermore, from Assumption 12.41 
for some positive constant C only depending on M, we have 

(47) T y {a*C,&*0) < CTV*O 1/2 ^(o-*0) 1/2 - 

Equation (l24l then reads, for any function ip € C, 
(48) 

Vx G R d , Ty(£(;y)x,&*D<p) = -(l/2)([o + H}(; y)x, Dip) 2 = -Ty(&*x,&*D<p). 
From (llOal ). (l46l ) and (l48l . we have for any function ip G C 

(x,A(y)x) = 2 lim T y (a*x + V°u x (-,y)x) = 2T y (a*x + £(■, y)x) 

A — ^0 

= 2T y {a*x + 1{-,y)x, a*x + a* Dip) 

< 2CT y {a*x + £{-,y)x) 1/2 T y (a*x + S*£><^) 1/2 . 

Gathering this with the inequality T y (a*x + a* Dip) < MM[|<x*a; + a*D<p\ 2 ] and d451) . we 
deduce < 2C 2 M(x,2x). 

It just remains to prove that the drift term B is orthogonal to K = KerA(y). Due to 
dlQcl ) and the fact that A' = Ker A(y) does not depend on y G M d , it suffices to prove that 
Ker H(y) C Ker^4(y) = K. But this is an easy consequence of (llOal) . (llObl) and Assumption 
1241 especially |ii"(a;, y)| < M 2 a(u, y). □ 
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10 Tightness 



We now turn to the tightness of the process X £ , ie we want to prove that the family (X £ ) £ is tight 
in C([0, T], M. d ) equipped with the uniform topology. That step of our result deeply differs from 
the uniform elliptic case Ifl3l . Indeed, uniform ellipticity of the diffusion matrix provides strong 
transition density estimates of the process X £ , the so-called Aronson estimates, from which the 
tightness of X s is then easily derived. Of course, in the degenerate framework, tightness of X e 
cannot be tackled this way. The method presented below is inspired from |[T5l and is based on 
the idea that the process X s is not too far from being reversible at a microscopic scale. The 
contributions of the macroscopic variations make a drift appear, unlike in IPT51 . 

Let us now go into details. As in Section [6l we can solve the following equation for i = 
1, . . . ,d and A > 

(49) Xw{(.,y)-S y w i x (.,y) = b i (.,y) 
and get the same properties as in Proposition [63j namely 

Proposition 10.1. For each fixed y G M rf and 1 < i < d, the famUy (V a w l x (.,y))\ converges 
to a limit Ci(-iU) G L 2 (Q) d as X goes to 0. The same property holds for the derivatives, 
that is, the families {V a d yj w x )\, (V 7 dy_. yk w x )\ (1 < i,j, k < d) respectively converge to 

dy 3 Ci(-,y)> ^ J% fcCi(-)y) L 2 (£l) d . Furthermore, the function w l x as well as its derivatives 
d yj w\, dy. yk w\ satisfy (16.31 ) and estimates (1221 ) and (1231) . for some positive constant C{]Q j\ 
independent of A > and y £ M. d . 

As in the proof of Theorem 18.21 we want to use a time reversal argument. Once again, we 
are faced with the lack of smoothness of w\ in order to apply the Ito formula. To overcome 
this difficulty, we proceed as in Section [7] Since the arguments are quite similar, we just outline 
the main ideas without further details. Let us consider, for n > 1, A > and 1 < i < d, the 
solution w x ,n of the following equation 

(50) Xw? (., y) - (, y) - n' l ^w\ n {.,y) = b t (., y) 

Introducing a sequence of regularizing sequence of mollifiers (ft n ) m eiN G C£°(M. d x 
(smooth functions with compact support), we define 

w \ m( w ' y)= w T( T x u > y - y')9m(x', y') dx' dy', 

which is a smooth function. Following the proof of Theorem 18.21 under the invariant measure 
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e 2V ( X ) dx of the process X n,£ , we can write 

(51) <' 2 yxr,X t "' £ ) =w^ n jX n /,K' £ ) + j\c^(w^J./e,.))](xr,X^)dr 



+ (Mf n ' m - M £ > n > m ) 



(52) w^jX n /,X^) =wfjX^' £ ,X^) + /Vt(«£J7^))]C^,^)*- 

Js 

+ (M £ t > n ' m -M £ < n > m ), 

where J\4 £ ' n,m and .A/f e,?1 ' m are two martingales respectively with respect to the forward filtra- 
tion (J r s' £ )o<s<T = <J {Xr ,e ; < r < s} and with respect to the backward filtration (G 1 s' £ )o<s<T 
a {Xr ,£ ; s < r < T}. The quadratic variations of both martingales match 

e ~ 2 f \Dw$ m + ed y w^ m T(a + n^ld)(Dw^ m + ed y w^J(T r ' £ , X^ £ ) dr. 
Jo ' 

Adding up (IBTI ) and (l52l . passing to the limit as m — ► oo (as explained in |[T3l Lemma 5.3]) and 
n — ► oo (as explained in Section Q and using d49l ) leads to 

(53) e" 1 / bi(X £ r ,X £ .)dr = f [ew^ + (l/2)tra,ce(aDd y wi 2 ))(X E r ,X £ ) dr 

J S J S 

rt 2V _ 

+ y —[divy(e- 2V a[Dw l e2 +ed y w l e2 ])]{X £ T ,X £ )dr 
+ (1/2) T Div(a) • d y w l £2 (X £ r , X £ ) dr 

J s 

+ e(M £ t -M £ s )+e(M £ t -M £ s ) 

= &S,t + At,t + I S,t + 1 S,t > 

where e.M e and e.M e are two martingales, respectively with respect to the forward filtration 
(J r£ )o<s<T = &{X £ ;0 < r < s} and with respect to the backward filtration (G £ )o<s<T = 
a {X s ; s < r < T}, with quadratic variations 

(54) I (Dw^ + edyW^faiDw^ + edyW^QC^ X £ ) dr. 
Jo 

Theorem [8j]establishes the following convergence 

limE £ [ sup El'l + El'l - [ G{X £ .)dr 

0<t<T ' ' Jo 
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where 

G(y) = M[(l/2)ti a ce(ad y i l )(.,y) + (e 2V /2)dW y (e- 2V a^)(.,y))]. 

From Proposition ! 10. H and (|22l . G is bounded so that the tightness of the process t \-t G(X £ ) dr 
in C([0, T], M) results from the Kolmogorov criterion. The tightness of E 1>e + E 2,£ follows. 

Let us investigate now the term T 8 f = (1/2) f s Div(a) -dyW 1 ^ {X 6 r , X £ ) dr. Note that it can 
not be treated with Theorem 18. II because the L 2 -norm of Div(a)d y w £ 2 need not be bounded. 
Inspired by the proof of Theorem [3J] in Section |9l we define 

= Div(a)-dywl 2 +tT&ce(aDdyW % e 2)+£diVy(a)-dywl2+etTace(adyyWl2)— 2eapjd yj u l £ 2d Vp V. 

By making two successive integrations by parts as in Section |9l we establish for any ip G 

C x C<f (R d ): 

PropMI „ . .„ 

My [** , cp] = -M v [ad y w l £2 • (Dip + edycp)] < C £ M V [\<r*(D<p + ed y ip)\ 2 } 1/2 , 
where the family (eC £ ) £ converges to as e goes to 0. Theorem 1 8 . 2 1 then ensures that 



sup ( / ^ £ {X £ r ,X £ r )dr 

0<s<t v Js 







-0<s<t 

as e goes to 0. Thanks to Theorem 18. H and Proposition llO.il we have 



sup | / trace(aDd y w £ 2 ) (X e r , X £ ) dr — / <&(X £ )dr\ 
o<s<t Jo Jo 



as e goes to 0, where $(y) = lim £ _ > oM[trace(ai^9yiM £ 2)(-, y)]. The Kolmogorov criterion 
and Proposition IIP. II ensure the tightness in C([0, t];M) of the process J Q $(X r £ )dr. More- 
over, from Proposition 1 1 0. 1 1 and (l30l) . the process f Q \ediv y (a) ■ d y w £2 + etr&ce(ad yy w £2 ) — 
2ea p jd y u l E 2dy p V] (X £ r , X s ) dr converges in law in C([0, T]; M) to 0. This proves the tightness 
of T 1 ' e 'inC([0,t];m d ). 

It just remains to treat the martingale term T 2,£ . According to Theorem 4.13 in (61, it 
suffices to establish the tightness of the brackets of these two martingales (see (l54l). Their 
tightness results from Theorem 18. II Proposition 110.11 and the Kolmogorov criterion again. The 
tightness of X s is now clear. □ 
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